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We present the current status of our computation of quark bilinear renormalization constants for 
Wilson fermions and Symanzik improved gauge action. Computations are performed in Numeri- 
cal Stochastic Perturbation Theory. Volumes range from 10 4 to 32 4 . 

Renormalization conditions are those of the RI'-MOM scheme, imposed at different values of the 
physical scale. Having measurements available at several momenta, irrelevant effects are taken 
into account by means of hypercubic symmetric Taylor expansions. Finite volumes effects are 
assessed repeating the computations at different lattice sizes. In this way we can extrapolate our 
results to the continuum limit, in infinite volume. 
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1. Motivations 

For a long time, lattice perturbation theory was the only available tool for the computation of 
Lattice QCD Renormalization Constants (RC's). Since the introduction of methods that allow the 
non perturbative computation of the RC's for generic composite operators [pig], these techniques 
are preferred. Nevertheless there is no theoretical obstacle to the perturbative computation of either 
finite or logarithmically divergent RC's. In principle, Lattice Perturbation Theory (LPT) provides 
the connection between lattice simulations and continuum perturbative QCD, that works only at 
high energy. The main difficulties of LPT are actually practical. First of all, LPT requires much 
more effort than in the continuum. Because of this, a lot of perturbative results are known only at 
one loop. Second, LPT series present bad convergence proprieties, so that we should emphasize 
that a fortiori one loop computations cannot give the correct results. A good example is provided 
by the computation of the quark mass renormalization constant: there appeal - to be discrepancies in 
between determinations coming from perturbative and non-perturbative techniques [||]. The very 
point is that as long as a comparison is made taking into account one loop, it is virtually impossible 
to assess the systematics of both results. 

In order to obtain higher orders in PT expansions, one can make use of Numerical Stochastic Per- 
turbation Theory (NSPT) [Q], a numerical implementation of Stochastic PT [Qj. Renormalization 
constants can be computed in NSPT at 3 (or even 4) loops, as it has been done in the case of Wilson 
fermions with non-improved gauge actions [Qj. A nice aspect is that one can work in the massless 
limit, where RC's are often defined, and no chiral extrapolation is needed. We will discuss how one 
can assess to finite a effects by means of hypercubic Taylor expansion. Moreover, finite volume 
effects can be corrected by repeating the computations at different lattice sizes. Measuring RC's 
for different values of rif is also possible to know the dependence on the number of flavors. 

2. RI'-MOM scheme 

The scheme we will adhere to is the so called RI'-MOM scheme, which became much popular 
since the development of non-perturbative renormalization. RI emphasizes that the scheme is reg- 
ulator independent, which makes the lattice a viable regulator, while the prime signals a choice of 
renormalization conditions which is slightly different from the original one. An important feature 
of this scheme is the fact that the relevant anomalous dimensions are known up to three loops in 
the literature [Q]. 

The main observables in our computations are the quark bilinears between states at fixed off shell 
momentum p: 

j dx{p\y(x)Ty{x)\p) = Gr{pa) (2.1) 

where Y stands for any of the 16 Dirac matrices, returning the S, V, P, A, T currents. Our notation 
points out the pa lattice space dependence. 

Since these quantities are gauge dependent we need to fix the gauge. We will work in Landau 
gauge, which is easy to fix on the lattice [j8|]. This gauge condition also gives an advantage: the 
anomalous dimension for the quark field is zero at one loop. 
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Given the quark propagator S(pa), one can obtain the amputated functions 

Gr{pa) — >Ty(po) = S~ l (pa)Gr(pa)S~ l (pa). (2.2) 

The Fr(pa) are then projected on the tree-level structure by means of a suitable projector Po r 

T {pa)=Tr(Po T T T {pa)). 

We can finally express the renormalization conditions in terms of the Or(pa) operator: 

Z 0r {na,g(a))Z-\na,g(a))O r {pa)\ p2=ti2 = 1. 

The Z's depend on the scale ji via the dimensionless parameter \ia, while the dependence on the 
coupling g(a) is given by the perturbative expansion. The quark field renormalization constant Z q 
in the formula above is defined by 



1 Tr^CH) 

A? — i 



' q ~ l 12 p 2 



(2.3) 



To obtain a mass independent scheme, we impose these conditions at the massless point. In the 
case of Wilson fermions this requires the knowledge of the critical mass. While one and two loop 
results are known from the literature (and can be used as a consistency test) the third loop is a 
byproduct of these computations. 

3. Finite lattice size effects 

At the generic n^Moop the RC's take the form 

n 

Z n = c» + dj(y)\og ] (pa) +F(pa), 

7=1 

where logarithmic contributions are known from the literature and one is mainly interested in the 
finite number c n . The first thing to do is then to subtract the divergent logs (again, we take them 
from the literature). After such a subtraction we still have the irrelevant contribution F(pa), that 
can be fitted by means of a hypercubic Taylor expansion. We show how this technique works by 
an example. 

Consider the two points function in the continuum limit: 

r 2 (p 2 )=s-\ P 2 ). 

On the lattice it depends on the dimensionless quantity p = pa. Furthermore, we explicit write the 
dependence on the coupling 

r 2 (p,m cn p- 1 ) = S- l (p,m cr ,fi- v ) 

= $ + m w (p) - t(p,m cr ,f}~ 1 ), 

where the hat stands for a dimensionless quantity. Here is the mass term generated by the 
Wilson prescription (it is G(p 2 )), t,(p,m cr , j3 _1 ) is the self energy and m cr is the quark critical 
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mass. The fermion mass counterterm arises because Wilson regularization breaks chiral symmetry. 
Both these last two terms are <^(j6 _1 ). 
The self energy can be written as 

t(p,m cr ,f5~ l ) = t c (p,m a -,f5~ l )+t v (p,m cr ,p~ 1 )+t other (p,m cr ,p~ 1 ). 

These are the different contribution along the different elements of the Dirac base. In particular, £ c 
is the contribution along the identity, Ly the contribution along the gamma matrices. 
The self energy E c contains the contribution of the mass: 

I(0,m cr ,j8 _1 ) = t c (0,m cr ,P~ ) = m cr = am cr . 

To show how hypercubic Taylor expansions work, we can concentrate on the £y term, from which 
we can in this way extract the quark field RC. We can consider a Taylor expansion of £y in the 
form 

t v (p,m cr ,r 1 ) = iZr^ (4 0) (A« C r,i3- 1 )+^4 1) (P>«-^- 1 )+^ 2) (A« C r,y3- 1 ) + . 
n 

(3.1) 

where the expansion entails an expansion in powers of a. The functions E^(/5,/n cr ,j3 _1 ) are com- 
binations of hypercubic invariants. As an example, the first term in the expansion ( |3.l| ) can be 
written as 

tPfaAcr,?- 1 ) = «i (0) l + 4° 5 I>v + afE^v + «f {LpI) 2 + M 

v v V v / 

As a general recipe, all the possible covariant polynomials can be found via a character's projection 
of the polynomial representation of the hypercubic group onto the defining representation of the 

group. One can see that plugging this expansion in the definition (23) of Z q the only term that 

fo) * (o) 

doesn't vanish in the continuum limit is G£j , the first coefficient of the expansion ofL v . 



4. Finite volume effects 

In the previous section we have shown how one can overcome the effects due to the discretiza- 
tion. These are not the only systematic effects in a lattice simulation. Though large lattices are 
available, one can't neglect the finite volume effects. This means that on a finite lattice we have to 
consider also a pL dependence in our quantities. 

We will now show how the analysis of £ = t.(pa, pL,am cr , jS -1 ) can be modified to take care of 
the finite volume effects. Let us take iLy\pa,pL,am cr ,fi~ l ) as an example. In the spirit of [0], 
consider the ansatz 

±®(pa,pL) = 4 0) M+ (t^(pa, P L)-t^(pa)) 
= £ v d \pa)+A±$\pa,pL) 

ty (pa, oo) and to keep the formula simpler we omit the dependence on m cr and 



where Ey (pa) 



4 



High loop renormalization constants for Wilson fermions/Symanzik improved gauge action M. Brambilla 



We can now expand as in (3.2): 



f$\p,pL) = 1 + of ZpI + «f 14 + «i 0) (ZpI) 2 + *^\PL) + (4-D 

V V V V / 

where the continuum limit ALy \pL) = ALy\pa = 0,pL) has been taken into account. The ratio- 
nale for this is the following: pL effects are present also in the continuum limit, and pa corrections 
on top of those are regarded as corrections on corrections. In principle, also these corrections can 
be taken into account, but this requires a larger number of parameters. 

The latter assumption has a strong consequence: measurements on lattices of different size are 
affected by the same pL effect once one consider the same tuples (n\ ,n2,ns,n4.) 

Innn 

PpL= —j—L = 2Tln yi . 

The practical implementation is the following. First of all one has to select a collection of lattice 
sizes and an interval \(pa) 2 min , (pa) 2 mx \- Then one considers all the tuples n = (n\ ,«2,«3,«4) which, 
for all the lattice sizes, fall in [(pa) 2 njn , (pa)l MX \ and chooses one representative for all the tuples 
connected by an H4 transformation. We further add to these data the measure taken at the highest 
value of (pa) 2 which falls in the interval \(pa) 2 min , (pa)l MX \ on the lattice with the biggest value of 
N = L/a. Assuming that this tuple (which we will call n*) is a good approximation to pL = oo, the 
measure will be considered as a normalization point. Is then possible to fit the parameters in the 
expansion (in our case the ct- ^), and the finite volume corrections, which came in the same number 
of the considered tuples. 

To be explicit, ( ft~l| ) can then be formulated in terms of the (n,N) dependence: 

tf\n ± n*,N) = a[ 0) l + oi 0) p 2 (n,N) + af/ (H,A0 + {p 2 (n,N)) 2 + Af + ff{a 6 ) (4.2) 

where n* means the tuple taken as normalization point. For the measure taken at the tuple n* the 
formula is slightly different: 

t { y\n*,N) = a[ 0) l + af ] p 2 {n*,N) + af ] p\n* ,N) + af ] (p 2 (n*,N)f + 0(a 6 ) (4.3) 



In (4.2) and (4.3) p 2 = Y*vPv an ^ P 4 = HvPv- Two issues of the procedure can only be assessed 
a posteriori: the assumption that the measure at the normalization point is free from finite volume 
effects and the stability of the fit. 



5. Results 

In the previous sections, we discussed the procedure to obtain results and keep finite a and 
pL effects under control. Since the work is still in progress, we only display the first technique 



at work. We computed at every order relevant expectation values given by (2.1) and amputated 
to obtain ( |2.2[ ). Finally we projected on the tree level structure and performed the pa analysis by 
means of the hypercubic Taylor expansion. 

We are not presenting any new result (apart a preliminar value for critical mass at three loop), but 
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we can compare some numerical results at leading order with analytical results [|H]] . 
The first quantity we consider is the simpler vector one: the measure of the RC for the quark 
propagator. This is a log free quantity, so we only need to extract the correct a — > limit. In 
figure |] we show the computation of Z q . It's easy to recognize the effect of different lengths of 
in the relevant direction. 

As an example of a logarithmically divergent quantity, in figure ||| we show the computation of 
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-0.58 - 




Figure 1: Computation at one loop of Z q . Different "families" of point correspond to the different compo- 
nents of momentum p along the relevant direction. In red measured values at given (p 2 ,Pn), in red the fit 
results and in blue the extrapolation at a = 0. 

Z s . In the figure, red point are the measurements before the log subtraction, blue the result after 
the y s log p subtraction. Since this is a scalar quantity, we haven't the effect of the length of the 
components in a given direction. 

In figure || we present the preliminar result for the critical mass at three loop, m$ . This quantity 
doesn't present either vector structure or log divergences. Critical mass is a byproduct of all the 
previous computations. The introduction of this counterterm is needed in the computation of the 
next loop quantities. 



6. Work in progress 

The aim of this work is to compute high order renormalization coefficients for masses, fields 
and bilinears for different actions. Analytic results are known at most up to 2 loops, while there are 
NSPT results for Wilson gauge-Wilson fermions up to 3 loops at various rif. We are in the process 
of taking into account pL effects for the latter action. We are now applying the method also to Tree 
Level Symanzik (current work) and Iwasaki gauge actions. For these actions we're also interested 
in computing the n f dependence. Current data are obtained from 32 4 lattices measured on an APE 
machine, while a C++ code is working on smaller lattices on standard workstations. 
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Figure 2: Computation at one loop of Z s . The anomalous dimension at one loop is different from zero, and 
the measure shows a logarithmic divergence (red points). In blue the results after the y s logp 2 contribution 
subtraction. 
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Figure 3: Preliminar result for the critical mass at three loop. In red the result of measurements, in black 
the result of the fitting procedure. 
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